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Despite being all roughly of solar composition, primitive meteorites (chondrites) present a diversity in their chemical, isotopic and 
petrographic properties, and in particular a first-order dichotomy between carbonaceous and non-carbonaceous chondrites. We 
investigate here analytically the dynamics of their components (chondrules, refractory inclusions, metal/sulfide and matrix grains) 
in protoplanetary disks prior to their incorporation in chondrite parent bodies. We find the dynamics of the solids, subject to gas 
drag, to be essentially controlled by the "gas-solid decoupling parameter" S = St/a, the ratio of the dimensionless stopping time to 
the turbulence parameter. The decoupling of the solid particles relative to the gas is significant when S exceeds unity. S is expected 
to increase with time and heliocentric distance. On the basis of (i) abundance of refractory inclusions (ii) proportion of matrix 
(iii) lithophile element abundances and (iv) oxygen isotopic composition of chondrules, we propose that non-matrix chondritic 
components had S < 1 when carbonaceous chondrites accreted and S > 1 when the other chondrites accreted. This suggests that 
accretion of carbonaceous chondrites predated on average that of the other chondrites and that refractory inclusions are genetically 
related to their host carbonaceous chondrites. 
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1. Introduction 

Oh 
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Q Meteorites are the witnesses of the first Myr of the solar sys- 
,H tern. Among these, chondrites have undergone little change 
^3 since their parent body (asteroid or comet) accreted, and thus 
i i preserve various solids formed in the protoplanetary disk: Most 
abundant in chondrites are chondrules, millimeter-sized silicate 
spherules that formed through rapid melting and cooling of pre- 
cursor material, via an as yet elusive mechanism ( |Connolly and| 
|Desch|2004 1. Metal and sulfide grains are probably genetically 
£f) related to chondrules (Cam pbell et a/.|2005| l. Refractory inclu- 
sions, in particular Calcium-Aluminum-rich Inclusions (CAI), 
■^j- likely formed by gas-solid condensation at high temperatures 
CD ( |Gros sman 2010 1, probably in a single reservoir located close to 
CN the Sun (MacPherson|2005|l — the "CAI factory" of |Cuzzi et~al. 
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All these components are set in a matrix of micron- 
sized grains. Radiometric dating suggests that CAIs are the 
'k^j oldest solar system solids, with an age of ~4568 Myr, while 
chondrule formation occured 1-3 Myr later (e.g. Bouvier anc 
C$ |Wadhwa|2010l|Villeneuve et aZ.|2009||Amelin et a/.|2010| l. 

While the abundances of nonvolatile elements in all chon- 



drites roughly match that of the solar photosphere (Palme and 



Jones 2005[), chondrites do present a diversity in their chem- 



istry, the isotopic composition of several elements, most promi- 
nently oxygen (Yurimoto et al. 2008), and in the textures and 
abundances of their different components (Brearley and Jones 
1998). 14 discrete chondrite groups have been hitherto recog- 
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nized, each of which representing either a single parent body or 
an array of parent bodies from the same reservoir. 

Meteoriticists have basically envisioned two possible expla- 
nations for this diversity: The first one is that the reservoirs rep- 
resented by each chondrite group acquired different composi- 
tions for their condensed fraction as a result of differing thermal 
histories. Little mixing is assumed between these reservoirs. 
Since differences between chondrite groups would be largely 
traced to what did and did not condense out of the initially hot 
gas, this "regional" view has come to be known as the incom- 
plete condensation model (e.g. |Wasson and Chou|1974 Bland 
\et a/.|2003] l. The second possibility is that mixing was extensive 
in the solar nebula, and that chondrite diversity results from the 
incorporation of essentially the same components, but in dif- 
ferent proportions depending on the group. Then, the different 
components need not be cogenetic. The latter would schemati- 
cally comprise a high-temperature component (chondrules and 
CAIs) and a low-temperature component (matrix). This is the 
two-component model (e.g. Anders|1964 Zanda et 01) 2006 ). 

Along with the mechanism(s) that caused fractionation — i.e. 
separation of compositionally distinct constituents — in chon- 
drites, it remains to be understood whether their diversity pri- 
marily reflects variations in space (i.e., heliocentric distance) or 
time. Rubin (201 1 ) illustrates the former view, with chondrites 
roughly arranged in order of increasing oxidation state with in- 
creasing heliocentric distance, ending with carbonaceous chon- 
drites. |Cuzzi et al.\ ( |2003| l and |Chamber"s| ( |2006| l advocate the 
latter view, with carbonaceous chondrites being the oldest. 

The age range spanned by the different inclusions within sin- 
gle chondrites suggests that they spent up to several Myr as 
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free-floating objects in the accretion disk. Significant fraction- 
ation may thus have resulted from their transport dynamics be- 
tween the times of their formation and their accretion. In most 
disk models, these dynamics are governed by gas-grain inter- 
action ( |Cuzzi and Weidensch illing 2006). This is supported by 
petro graphic evidence: Chondrules and other non-matrix chon- 
drite components in a given meteorite exhibit a narrow distri- 
bution as a function of the product of the density and the ra- 



dius, which is the relevant aerodynamic parameter ( Cuzzi and] 
|WeidenscMiring|2006||Kuebler et aZ.p9l?9l|May etal.\\999) . 
Aerodynamic sorting could also account for metal- silicate frac- 
tionation (see e.g. Zanda et al. 2006). Finally, the small size 



(< 20 /vm) of high-temperature minerals identified in cometary 
material is consistent with the increasing difficulty of outward 
aerodynamic transport of inner solar system material with in- 
creasing size (Hughes and Armitage 2010]). 

Several studies have addressed the dynamics of solids within 
the disk (Cuzzi and Weidenschilling 2006), with focus on re- 
fractory inclusions or other high-temperature materials (e.g. 
IBockelee-Morvan et aZ.|2"002"l|Wehrstedt and Gail|2002t |Cuzzi| 
\et a/.|2003|[Ciesla|2009a|bl[20T0lpighes and Armitage |2010), 
the incomplete condensatio n model (|Cassen|1996||2001[ |Ciesla 
2008 ), water transport (e.g. Ciesla and Cuzzi|2006| )~ 



In this paper we revisit the radial dynamics of chondrite com- 
ponents from a simplified, analytical point of view, taking the 
above works as guides. Our goal is to find out whether the di- 
versity of chondrites can be understood in terms of the transport 
of their components in the disk, adopting standard modeling as- 
sumptions. We will show that inferences can be made in this 
respect independently of the details of the structure or turbulent 
properties of the disk, which are still uncertain. 

The outline of this paper is as follows: In Section [2] we 
briefly review a first-order classification of chondrites in car- 
bonaceous and non-carbonaceous chondrites. We review the 
dynamics of gas and solids in Sections [3] and [4] respectively. 
In Section [5] we introduce the "gas-solid decoupling param- 
eter" S central to this study and justify an approximation of 
the one-dimensional, vertically integrated continuity equation. 
In Section [6] we present four arguments supporting the con- 
jecture that non-matrix chondrite components had S < 1 
when carbonaceous chondrites accreted and S > 1 when non- 
carbonaceous chondrites accreted. Cosmochemical implica- 
tions are discussed in Section|7] In Section|8] we conclude. 




Figure 1 : Photomicrographs of thin sections of a carbonaceous chondrite (Al- 
lende, top) and an ordinary chondrite chondrite (Hallingeberg, bottom), viewed 
in polarized light. The round objects (sometimes fragmented), generally with 
colorful crystals, are chondrules. Irregularly-shaped, fine-grained objects with 
bluish-whitish shades visible in the upper photograph (mostly near the center) 
are CAIs. All these objects are set in a matrix (black in transmitted light). 



2. A first-order classification of chondrites 

It is certainly beyond our scope to propose particular geneses 
for each chondrite group. Instead, we shall content ourselves 
with a first-order classification of chondrites in two super-clans 
( Kalleme y n et al. | 1 9 96 ; Warren 2011 1: the carbonaceous chon- 
drites and the non-carbonaceous chondrites (see Figure [T] and 
Table [T). 

The carbonaceous chondrites (hereafter the "CCs") com- 
prise the CI, CM, CO, CV, CK and CR groups— we leave 
aside the CH and CB chondrites, which may have formed in 

Krot et al.\ ( |2005| ) but see 



very peculiar environments (e.g. 



Gounelle et q/.|[2007 ) — and are the most chemically primitive 



Table 1 : Compared primary properties of chondrites 



Properties 



C chondrites EOR chondrites 



CAIs 


abundant 


rare 


matrix fraction 


high 


low 


monotonic trend of lithophile 






elements with volatility? 


yes 


no 


Mg/Si 


solar 


subsolar 


chondrule O isotopes 


ls O-rich 


16 0-poor 




variable 


less variable 
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chondrites. The non-carbonaceous chondrites comprise the 
£nstatite (EH, EL), the Ordinary (H, L, LL) and the /?umuruti- 
type (R) chondrites and will be henceforth be referred to as 
"EORs". Schematically, the main distinguishing features are 
(see |Brearley and Jones| ( [1998) ; |Scott and Krot| pOOl) for de- 
tailed descriptions of chondrites): 

(i) CCs have higher CAI abundances (up to > 5 % in volume; 
e.g. |Hezel et af.|2008] ) than EORs (<0. 1 %). 

(ii) CCs have a higher proportion of (fine-grained) matrix 
(generally >30 %) than EORs (> 10 %). Chondrule textures 
also testify to higher dust concentrations in their formation en- 
vironment for CCs (Rubi np"010[ ). 

(iii) CCs exhibit a depletion (compared to solar abundances) 
of moderately volatile elements relative to refractory elements, 
as an essentially monotonic function of the 50 % condensation 
temperature, while EORs exhibit no such clear correlation (ex- 
cept for the most volatile elements). 

(iv) The Mg/Si ratio is solar for CCs and subsolar for EORs. 

(v) CCs have more 16 0-rich chondrules than EORs. Also, 
in the three-isotope diagram, chondrules from single CCs de- 
fine larger fields, which partly overlap among different chemi- 
cal groups, while chondrules from single EORs form narrower, 
resolved fields ( |Rubin|2000HClayton|2003||Jones et al.\2000) . 

Regarding CCs, while properties (iii) and (iv) are true of 
them as whole rocks, they do not seem to hold for their indi- 
vidual components sep arately (|Bland et a/.| ( [2005] l; Hezel and 
Palme] ( |2010| ), but see |Zanda et al. |2012 1. This is known as 
complementarity. 

A primary assumption in this work will be that that the pro- 
portions of components in chondrites reflect their relative abun- 
dances around the disk midplane prior to accretion. It is how- 
ever conceivable (e.g. Ciesla 2010) that the accretion process 
(if not gravity-controlled) caused preferential incorporation of 
components of certain sizes: This we shall refer to as the "ac- 
cretion bias". A prototypical example is the turbulent concen- 



tration scenario of Cuzzi et al. ( 2001 1, but similar sorting calcu- 
lations would be desirable e.g. for the streaming instability or 
for "hit-and-stick" growth. The evidence for complementarity 
in CCs would militate against any accretion bias between chon- 
drules and matrix grains. At any rate, the lack of a systematic 
chondrule size difference between CCs and EORs suggests that 
accretion bias would not account for the differences between 
the two super-clans as wholes. 



3. Disk modeling 

The disk is described in a cylindrical coordinate system, with 
R the heliocentric distance, z the altitude above the midplane, 
and the azimuthal angle. We denote by u, p, T and P = pc 2 s 
the gas velocity, density, temperature and pressure, respectively, 
with c s - y/ksT/m the isothermal sound speed, where kg is the 
Boltzmann constant and m the mean molecular mass (here 2.33 
times the proton mass). £2 and v K = OR are the Keplerian 
angular and linear velocities, respectively. 

We treat the disk as vertically isothermal (see Appendix A), 



and thus vertical hydrostatic equilibrium implies: 



p(R,z) = 



exp - 



(1) 



with the pressure scale height H = c s /Q and the surface density 

We assume the disk to be turbulent and assume axisymmetry 
in the sense that variations of any quantity Q in the azimuthal 
direction may be treated as turbulent fluctuations about a mean 
Q. We denote the Eulerian perturbations with 6Q = Q — Q. 

In that case, from angular momentum conservation: 
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1 d IR 112 , 

R^dR\~cr j 
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VK 



(2) 



with T R( p and T n 



the Rip and zcf> com ponents of the turbulent 

and a a dimen- 



stress tensor (Balbus and Papaloizou 



1999 



sionless number parameterizing the former as: 
3 



I Rip 



-aP 



(3) 



It must be noted that defined as such, a may a priori vary in 
space and time, and in particular, we do not make the assump- 
tion that the stress tensor obeys a viscous prescription (e.g. 
Takeuchi and Lin|2002 1. Observation s of disks suggest a v al- 



ues around 10 2 at 10-100 AU scales dHartmann et al. 1998 1, in 



the range of those achieved by the magneto-rotational instabil- 



ityjBalbus an d Hawley|1998||Lesur and Longaretti|2007[ |Bai 
|and Stone|201 1 1. In the dead zone ( Gammie|1996| l, where the 
latter is unlikely to operate, a should drop to 10 6 - 10 A around 
the midplane (e.g. Fleming and Stone[2003 ; Ilgner and Nelson 



2008] |Oishi and Mac Low|2009||Turner et q/.|2010) . 
The mass accretion rate can then be expressed as: 



■I 



M = -2kR pu^dz = 6ttR 



1/2 



d_ 
BR 



where 



v 41 



c c 
a^pdz = {a) p 



Q 



(4) 



(5) 



with (a) p the density-weighted vertical average of a. 



Note that we hereby ignore infall (unlike e.g. |Zhu etal.\2010) 
and photoevaporation (unlike e.g. |Desch|2007[ [Hughes an d Ar-| 
mitage 2010), which may be most important in the early and 
late stages of the disk, respectively. However, while these ef- 
fects would influence the evolution of S, they would not affect 
the validity of the radial velocity formula |2]i in the disk interior. 

An important benchmark for this paper is the steady-state 
disk, where M is uniform. This should be a good approximation 
of the inner regions of a real disk if the local evolution timescale 



R 

t vis (R) = — =0.2Myrfl 



1 km/s\/10 



(6) 



'Noting B the magnetic field, they are defined as (i = R, z): 



T i( j, = pSujyduj 
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(which is the timescale of accretion to the Sun of gas at helio- 
centric distance R) is shorter than the disk's age. In that case, if 
neither S nor v diverges as R — > 0, we have: 



M = 3ttEv 



4. Dynamics of solid particles 



(7) 



We now proceed to the transport of solid particles. By "parti- 
cle", we mean any detached solid body, belonging to any of the 
chondrite components mentioned in Section [T] We first review 
the dynamics for a single particle, before writing the mean-field 
continuity equation for a population of identical particles. 

4.1. Single particle dynamics 

The dynamics of a small particle are dominated by gas drag. 
For a spherical solid particle of radius a and internal density p s , 
the stopping time is (Epstein 1924) 



n p s a 
8 pc s ' 



(8) 



For a 0.3 mm-radius chondrule, p s a = 1 kg/m , which we take 
as our normalizing value. Within a factor of a few, it indeed 
pertains to most chondrules, metal/sulfide grains and refractory 
inclusions alike in most chondrite groups (King and King 1978 ; 
Brearley and Jones|1998||Kuebler et aZ.|1999||May et a/.|1999| l. 



We do not study here the fate of larger (> 1 cm) bodies. This 
is because (unbrecciated) chondrites are homogeneous at the 
centimeter scale and thus their composition was determined at 
the agglomeration of their millimeter-sized component^] 

A measure of the coupling of the particles to the gas on an 
orbital timescale is the dimensionless stopping time 



St = Qr 



n p s a 
2~1T 



2x 10" 



p s a \/10 3 g/cm 2 



1 kg/m 2 



(9) 



where the second equality holds at the disk midplane. Clearly, 
St <K 1 for millimeter- or centimeter-sized particles and smaller. 
In this limit, gas-solid drag forces the particles to essentially 
follow the gas, but with a systematic relative drift due to the 
pressure gradient, so that the particle velocity is given by 
( |Youdin and Go odman 2005): 



VP 

P ' 



(10) 



The radial component of the second term on the right-hand side, 
rdP c] 

Vdrift = -"5= ~ -St , (11) 

p oR v K 

is robustly negative in disk models (except perhaps around the 
inner edge of the dead zone, e.g. Dzyurkevich et a/.|2010| , and 
alone entails a sunward drift on a timescale 



^draa — 



R 



|Vdri 



OSt I c 5 



= 2My<£ 



1 km/sr/ 10- 



St 



4.2. Averaged continuity equation for a fluid of particles 

We now consider a population of identical, non-interacting 
particles and treat it as a fluid of density p p (not to be confused 
with the internal density p s of each particle). The continuity 
equation averaged over turbulent fluctuations may be written in 
the form ( |Schrapler and Henning|2004}|Balbus|201 1) : 



dp P 
dt 



ld_ 
RdR 



R P P v p ,r - DrrP 



dR 



m 



d_ 
~dz 



p p v p , z - Dzzp 



d_(pp 
dz\p 



= 0. 



(13) 



We shall henceforth drop the overbars. For St <k 1, the diffu- 
sion coefficients equal those of passively advected scalars in the 
gas ( Youdin and Lithwick|2007 1 and will be parameterized as: 



Drr - 8 R - 



and D 72 = 5 7 



(14) 



Q Q 
For correlation times of order Qr l (Fromang and Papaloizou| 



2006; Cuzzi and Weidenschilling 20061, 6r and 5 Z should be 



comparable to a because of the coupling between azimuthal and 
radial components due to Coriolis forces. We coin the radial 
and vertical Schmidt numbers^ 



Sc s = 



and 



Sc z = 



(15) 



|Prinn| ([T990jl and |Dubrulle and Frisch| ( fTWT) argued that, in 
hydrodynamical turbulence, Scr should not exceed 1, as is the 
0.176 value obtained in turbulent rotating flow experiments by 



|Lathrop et al.\§\992\ or the ~ 0.7 values mentioned by |Gail 
( |2001| l. |Johansen et al\ ( |2006| ) reported Sc s = 4.6a 26 and 



Sc, = 25.3ar U4b from local MHD net-flux simulations, which 
thus can exceed 1, consistent with most earlier results ( Joh ansen] 
\et a/.|2006) l. These uncertainties make it prudent to keep track 
of the Schmidt numbers throughout this study, although we gen- 
erally expect them to be of order unity. 

5. Net radial transport of solid particles 

We now wish to derive the vertically-integrated continuity 
equation of the particle fluid. Beforehand, we will introduce 
the dimensionless numbers of special relevance in this problem 
and briefly review the vertical distribution of the particles. 

5.1. The gas-solid decoupling parameter S 

The ratio of particle-gas drift and gas advection contributions 
in the radial transport is, using equations (|2]i and ( 1 1 

Vdrift St _ _ np s a 
ur a 2 Ea 

p s a \/10 3 \/10 3 g/cm 2 



(16) 



0.2 



.(12) 



lkg/m 2 /\ a 



2 However, some of these larger bodies may have (i) travelled to another 
reservoir and (ii) contributed by shattering to its sub-millimeter-size inventory. 
This is thus ignored here. 



There is some variation in terminology in the litterature: This number 
is called a "Prandtl number" by Prinn ( 19901 and Cuzzi and Weidenschilling 
1 2006i, and Youdin and Lithwick ( 2007 1 use "Schmidt number" for the ratio of 
the diffusivity of gas to that of the particles (1 for St <sc 1). 
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the last two equalities holding at the midplane. We also define: 

St St 
S R s — =SxSc s and S-= — = S x Sc 2 (17) 

or o z 

As we shall se e in the next su bsection, S z is a measure of verti- 
cal settling (Cuzzi et al. 1996 and we will see later in Section 
|6. 1 1 that Sr characterizes the ratio of the particle-gas drift and 
the diffusion contributions to the radial flux. 

So long the Schmidt numbers are of order unity, S , Sr and S,- 
actually are of the same order of magnitude. Hence, the qualita- 
tive aerodynamics of the particles are essentially controlled by 
one single number, namely S . Schematically, finite-size effects 
can be ignored for S « 1, while particles significantly decou- 
ple from the gas for S » 1 . In that sense, one may call S the 
"gas-solid decoupling parameter". To keep track of the exact 
role of S , Sr and S z , we shall however remain careful in writ- 
ing either one in the formulas (unless otherwise noted, they will 
be implicitly understood as being evaluated at the midplane). 



5.2. Vertical distribution of particles 

The vertical flux provides the dominating balance in equation 
( fl"3j ). On timescales longer than the vertical mixing timescale 



1 



Q max((5 z , St) ' 
the vertical distribution of the particles obeys: 

d . P P c . z 
T-ln— = -S z {R,z)—t, 
oz p H l 



(18) 



(19) 



from which the role of S z as a measure of settling is obvious, 
with dimensional analy sis indicating a dust layer thickness of 
order H/ yjl +S Z (e.g. 



Dubrulle et al. 



1995). 



5.3. One-dimensional continuity equation 

We now wish to integrate the two-dimensional continuity 



equation ( 1 3 1 over z to obtain the evolution equation of the par- 

X+co 
Ppdz. We must thus calculate 

the p p -weighted average of the particle velocity, similarly to 

Takeuchi and Lin ( 2002| l. For any quantity Q and any weight 

function w(z), we define its w- weighted average as: 



IT 6(zMz)dz 
£ + Jw(z)dz 



(20) 



Vdrift(z) can be readily expressed from equations (Hi and ([TJ 



Same does not hold true, however, for u R . As it stands, equation 
([2]i expresses it in terms of turbulent correlations whose depen- 



dence on z are essentially unknown (e.g. Miller and Stone 2000 
From ang et a/.||20lT i. Moreover, the nature of the turbulence 
may change with height (e.g. Fleming and Stone 2003] [Turner] 
et a/.| |20i0] l. In the current state of accretion disk theory, the 



vertical flow structure is thus uncertain. 



This is what Cuzzi et al. | jl996| called S (they had taken S' 



However, as we shall show, the net flow of the particles is 
actually weakly sensitive to the vertical gas flow structure, pro- 
vided Scr ~ Sc, ~ 1. In the limit S z » 1, the particles concen- 
trate around the midplane so that 



<v P ,R>p p ~ u R (R,0) + v drif ,OR,0) * Vdrift(*,0), 



(21) 



where the latter approximation uses the scaling Vdrin7«« ~ S 
(and hereby the hypothesis that Sc z ~ 1). In the limit S z <K 1, 
particles are well mixed vertically, and we have 



<Vp,R>p p ~ (u R ) Pp ~ {u R ) p = - 



m l i 2 dR 



[r^v). 



We thus see that the approximation 



<Vp,R>p p ~ (Ur) p + V dli ft(R, 0) 



(22) 



(23) 



should be accurate in either limit. For S ~ 1, the approxima- 
tions used break down, but, unless Sc z is very different of 1, 
or (a) p greatly exceeds a near the midplane, the imparted error 
should not exceed a factor of a few then. 

As regards the diffusion term, vertical integration yields: 



■r 



d f pdz 



(24) 



where / is defined by 

P P (R,z) W 
p(R,z) ~ W) 



f(R,z). 



(25) 



In principle, the second term on the right-hand-side of equation 
( 24 1 should be included as an effective correction to (v Pi r) p in 
the net flux. 

However, this correction is negligible compared to Vd r jf t (/?, 0) 
if S » 1 (as it is then of order (u R ) p ) and negligible compared 
to (u R ) p if S <K 1 (as / « 1 so that df/dR « 1). We shall thus 
drop it from the vertically-integrated flux. 

Hence, the one-dimensional continuity equation sought is: 



-^ + ^ M [R^ P ((uR) P + v <sm (RM 



-<Ak>p,s|j(^ 111 = 



(26) 



We will henceforth drop the {...) pfi and refer to Vdrift(^, 0) sim- 
ply as Vdnft- This equation is essentially that used in one- 



dimensional models in the literature (e.g. Cuzzi et al. 2003 



Ciesla 2010; Hughes and Armitage 2010J. We have shown here 
that it is accurate (within order-unity corrections, up to a factor 
of a few for S ~ 1), quite insensitively to the as yet uncertain 
details of turbulence modeling, and thus most appropriate to 
evaluate net radial transport. 
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Table 2: Properties of S 



Property 


S < 1 


S > 1 


Section 


Significant outward transport? 


yes 


no 




6.1 




Settling of mm-sized bodies? 


no 


yes 




6.2 




Coherence between grains? 


yes 


difficult 




6.3 




Radial diffusion 


unlimited 


limited 




6.4 




Proximity to Sun/epoch 


close/early 


far/late 




7.1 





6. The value of S at chondrite accretion 

In this section we shall advocate the following conjecture: 

Non-matrix chondrite components had S < 1 at the ac- 
cretion of CCs and S > 1 at the accretion ofEORs. 

Since S depends on size, we stress that we are speaking 
here of the largest, compositionally distinctive chondritic 
components, that is chondrules, CAIs, metal and sulfide grains. 
As mentioned in the introduction, the range of p s a of these 
inclusions in single chondrites is narrow, and roughly indepen- 
dent of the nature of the inclusions (e.g. King and King| 1978 



Kue bler ef"aT]|1999| |May et al\\999) , such that speaking of 
one value of 5 for a given chondrite should not be ambiguous. 
Moreover, mean chondrule sizes of different chondrite groups 
do not vary by more than a factor of a few (Bre arley and Jones] 
1998 ) so that, given the order-of -magnitude variations of L and 
a, we can even schematically think of all these inclusions in 
all chondrites as having the same p s a as 1 kg/m 2 . Therefore, 
variation in the value of S at the accretion of the different 
chondrites will essentially reflect differences in gas properties 
rather than intrinsic variations of the acccreted components. 

The above conjecture is supported by four arguments which 
we shall now examine each in turn. In each subsection, we 
begin with a property of S (all listed in Table [2]i and then link it 
to one of the meteoritic observations of Sectionp2j 

6.1. Outward transport and CAI abundance 

We wish to show here that particles formed close to the Sun 
cannot be transported much beyond the S — 1 line in significant 
amounts, be it through outward advection or outward diffusion. 

As regards advection, the argument essentially relies on the 
scaling relationship ( 17 1 Vdrift/ ur ~ S mentioned in Section 5.1 
For S » 1, the gas-drag induced drift velocity dominates and 
is robustly inward, and thus no significant outward transport 
is possible in this region, even if the net gas flow is outward 
(e.g. |Desch|2007 Jacquet et a/.|201 1 1. Thus, if a given region 
contains a significant abundance of a given component with S > 
1, the latter cannot have been transported from the vicinity of 
the Sun to there but must have formed in situ or further out, 
unless it was delivered by some mechanism other than the in- 
disk transport envisioned here (e.g. via jets as in Shu et t 
2001). One illustration are the "midplane flow" simulations of 



Hughes and Armitage (2010): for example, in their Fig. 5b2, 
20 yt/m particles pile up at ~90 AU, corresponding to S = 3.2. 



Turbulent diffusion would not alleviate this constraint. To 
s>ffset particle-gas drift, it indeed requires the particle-to-gas ra- 
.tio to vary on a lengthscale R/S r, amounting to a steep inward 



gradient if Sr » 1. For example, Hughes and Armitage ( 2012| ) 
found in their simulations that particles were confined to re- 
gions at most with St < 10 _1 , i.e. S < 10 for their a = 10~ 2 . 

It is instructive to calculate the particle-to-gas ratio profile 
_that results from balancing the different contributions of the net 



radial flux of the particles in equation ( 26 1, that is: 



d (X p \ 

(Mr + v<jrift)Sp - DrrZ— I — | = 0. 



(27) 



From this, we obtain: 



d_ 
dR' 



JL 



ln^ = 



Ur + Vdrift 



D 



RR 



= -3Sc R ^ln(/? 1 / 2 2v) + 5 R ^.(28) 



dR 



where we have injected equations (|4]) and ( pT| . Assuming Scr 
to be constant, this may be integrated to yield: 



«p(J*s*$F<w) 



(ZvR 1 / 2 ) 



3Sc* 



(29) 



If we take the surface density of the gas as a proxy for the total 
abundance of solids, Ep/E is proportional to the fraction of the 
component in question (say, CAIs) in chondrites at their accre- 
tion location. 

To be more definite, we specialize to the case of a steady- 
state disk and assume radial dependences of the form X oc R~ p 



and T oc R ?. Equation (29 1 becomes: 

f exp(-cS)^ 
~rW 



Sen 



(30) 



with c = (2p + q + 3)/(3 - 2q) a number of order unitjj^] This is 
plotted in Fig. For S <*c 1, the R- 3Sc «/ 2 dependence matches 



that found by Clarke and Pringle (1988) for a passive contam- 
inant. The new result here is that Sr = 1 clearly represents a 
cutoff for outward diffusion, but depending on Scr (see Section 
|4.2| >, diffusion of inner solar system material may not necessar- 
ily have been efficient until even there, especially if the helio- 
centric distance must climb one or two orders of magnitude. 

The S — 1 line thus schematically represents the maximum 
range of outward transport within the disk. It is set by the 
disk model and may be located a priori, before any simulation. 
Then, the high abundance of CAIs — which formed close to the 
Sun — in CCs relative to EORs implies S < 1 for the former. 

6.2. Settling and proportion of matrix 

We mentioned in Section [572] that S z was a measure of set- 
tling. While for S z <sc 1, the particles may be considered to 
be well-mixed vertically, for S z » 1, particles are concentrated 
around the midplane, with their particle-to-gas ratio being en- 
hanced by a factor of order -\/S~ z (Dubrulle et al. 1995 1. 



5 We have used equation j37[ . If a is constant (and thus p + q = 3/2), c = 
17 74 and 11/4 for the viscous-heating- and the irradiation-dominated regimes, 



respectively (see Section [TT| . 
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Diffusion in steady disks 
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Figure 2: Plot of the equilibrium diffusion profile of CAIs (or any chondrite 
component produced near the Sun) in a steady disk. The "normalized CAI 
fraction" is proportional to the CAI-to-gas ratio and is set to unity for T > 1500 
K. Results are shown for two mass accretion rates, 10~ 7 (red) and 10~ 8 (green) respectively), but are spatially resolved after a timescale f sep . If accretion has 



Heliocentric distance 

Figure 3: Sketch of the decoherence of a chondrule population (red) and a dust 
grain population (blue) as a result of differential drift. We plot their probability 
density as a function of heliocentric distance after a chondrule formation event 
at time t = (when these probability densities are essentially delta functions). 
Due to turbulent diffusion, the two populations largely overlap for small t de- 
spite a faster mean drift of chondrules (here velocities are denoted v cn and Vd u8t , 



M Q .yr . Solid lines assume Sen = 1 and dashed ones Scr = 1/3 (an enhanced 
radial diffusivity). Vertical dotted lines mark the heliocentric distance where 
5 = 1 for both values of M. We have taken fj = 1, K = 5 cm 2 /g and a = 10~ 3 
(see SectionjTTJ. 



Thus, the abundance of chondrules with S z » 1 is en- 
hanced relative to comparatively well-mixed micron-sized dust 
grains. The maximum enhancement is reached when the 
dust grains themselves start to settle, i.e., when S Zt d ust = 
(0°jfl)dust/(Psfl)ch)5 c > 1, where the "dust" and "ch" sub- 
scripts refer to dust grains and chondrules, respectively: as- 
suming similar densities for both, this maximum enhancement 
is V a ch/fldust, which, for micron-sized matrix grains and chon- 
drule radii of a few tenths of millimeter (see Section 4.1), eval- 
uates to ~10. In view of this, the fact that EORs tend to have 
less matrix (< 10 %) than CCs (> 30 %) would be accounted 
by 5 ; < 1 for CCs and S z » 1 for EORs. 

Before closing this subsection, we comment that S- > 1 for 
EORs would also be consistent with a lower absolute dust con- 
centration in the formation region of their chondrules. Indeed, 
for a fixed 6 Z , a larger level of settling would correspond to a 
smaller column density of the gas, and hence of the solids for a 
given dust-to-gas column density ratioj^] 

6.3. Coherence of a grain population and complementarity 

Let us consider two particle populations, labelled 1 and 2, 
with different sizes, both concentrated at the same given helio- 
centric distance R at t = 0. This is sketched in Fig. [3] Although 
the differing sizes imply different mean velocities, these two 



not taken place by then, any chondrule-matrix complementarity will be lost. 

populations do not instantaneously decouple from each other, 
because turbulent diffusion widens these distributions, such that 
they continue to overlap at short times. We call coherence be- 
tween two populations this state where both populations are not 
spatially resolved from one another. However, as the diffusional 
spread goes like yfi while the separation between the two peaks 
increases linearly with t, this coherence is expected to be lost 
by a timescale f sep where the peaks are mutually resolved, and 
which we define with 



kdrifU - Vdrifulfsep - 2 ^2D SR f sep , 

from which we deduce 

tsep = * DRR [( Tl -£)dP/dR 



(31) 



(32) 



If the size of population 1 particles is much bigger than that of 
population 2 particles, e.g. if the former are chondrules, or re- 
fractory inclusions, while the latter are matrix grains, we have: 

,1/2 



r sep 



St 2 



2Myr 

^drag 

~Sr~ 



Sr 

tv 



10-" 



1 km/s ' 



S 2 Sc R 



(33) 



6 This would not be compensated by concentration at the midplane if .S dust > 
1 because then, defining 6d ust = X aust /X, the dust density there is 



\^z,dust 



41nB 



gust I £ 



where St, S and S r pertain to the larger particles. 

For X « 1, f sep is longer than t vls and thus "decoherence" 
cannot occur within the accretion timescale (which is necessar- 
ily bounded by that of the drift to the Sun). In this case, chon- 
drules and dust grains remain coherent. Note that this is meant 
in a statistical sense, that is, the two populations as wholes are 
not spatially resolved from one another ; of course, a particular 
chondrule and a particular dust grain floating nearby would not, 
on average, remain together unless they quickly coagulate. 

In contrast, for Sr » 1, f sep is smaller than f dra g, and judg- 



which, everything else (e dust , H, d z ) being equal, decreases with decreasing E. ing from equation (33 i, decoherence then seems unavoidable 
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if chondrules and matrix do not agglomerate within <k1 Myr 
in an evolved disk. That the accretion timescale was longer is 
evidenced by the broad chondrule age distribution in single me- 



teorites (Villeneuv e et aZ.|2009 1. In the "inward transport" sim- 
ulations of Ciesla ( 2009b| l, where particles are initially seeded 
at 20-30 AU, significant decoherence between the a = 0.5 mm 
and the a — 5 fim particles can be observed for M - 10~ 8 M Q /yr 
(his Fig. 18), where the S = 1 line lies at R = 4 AU. 

Matrix-chondrule complementarity in CCs is evidence for 
coherence between chondrules and dust grains. Indeed, were 
chondrules unrelated to the matrix, that the rock resulting from 
their agglomeration is solar (e.g. in terms of the Mg/Si ra- 
tio) would seem too coincidental (e.g. Hez el and Palme|2010 1. 
Therefore, chondrules must have formed among dust grains 
now found in matrix, with exchange of chemical elements tak- 
ing place between them during the chondrule-forming event 
while preserving a solar bulk composition. In addition, for this 
complementarity to have been maintained in CCs, both popula- 
tions of objects must not have drifted apart until accretion: from 
the above discussion, this suggests that formation and accretion 
of CC components took place in regions where S < 1 . 



Following Cuzzi et al. (2005), we stress that the source of 
the components of a given CC needs not be unique: If two 
source reservoirs A and B are solar in bulk composition, their 
respective contributions to a chondrite-forming region C will 
be solar too — since transport is unbiased with respect to com- 
ponent size (for S <s. 1) — and so will be the chondrite resulting 
from their mixing in region C. This chondrite will thus exhibit 
chondrule/matrix complementarity. Such an "hybrid" comple- 



mentarity may be exemplified by CM chondrites, which Cassen 
( |1996[ l suggested to result from the mixing of CV/CO-like with 
Cl-like materials. Hence, complementarity in one chondrite can 
be the result of several distinct chondrule-forming events. 
As regards EORs, their composition was interpreted by 



Larimer and Anders ( 1970) as resulting from a loss of refrac- 



tory, olivine-rich components from a starting material of solar 
composition. Lost components may include amoeboid olivine 
aggregates — a variety of refractory inclusions — ( |Ruzicka et al. 
|2012| ) and chondrules, which tend to have supersolar Mg/Si in 
CCs ( |Hezel and Palme |2010| ). Such a fractionation requires 
S > 1, again consistent with our conjecture. 

6.4. Turbulent mixing and chondrule oxygen isotopes 

Let us now focus on a single population of particles concen- 
trated at some heliocentric distance at some initial time t — 0. 

For S <k 1, diffusion will smear the particle distribution ac- 
cross the whole range spanned in the inward transport. Indeed, 
using equation ([6|, the eventual spread is of order: 



R 



(34) 



|Clarke and PringTe] ([1988) derived an exact solution illus- 
trating this in the case of a steady disk with S oc R 2 
and constant Sc R , for perfectly coupled particles: The par- 
ticle distribution in terms of heliocentric distances is a log- 
normal, which, expressed as a function of InR, peaks at R = 



Ro exp (-3v(Ro)t/2R(i) (with Rq the initial heliocentric distance) 
and has standard deviation cn nR = 2 -^v(Ro)?/Scr/?q. For 
t ~ f v ; s , the latter is of order unity, implying smearing of the 
distribution over a lengthscale comparable to R by then. 

In contrast, for S » 1, drag-induced inward drift over- 
whelms diffusion. Using equation ( 12 1, we haveQ 



ra 



R 



rag 



(35) 



Thus, a reservoir with S » 1 will not have received sig- 
nificant contributions having diffused from further than this. 
Hence, for a given size, mixing between regions separated by 
more than R/ ~\/S~r is insignificant. 

This may explain why the oxygen isotopic composition in 
chondrules exhibits less variability in individual meteorites for 
EORs than for CCs, if we have S > 1 for EORs and S < 1 
for CCs. It is possible that the differences between the three 
ordinary chondrite groups — a subset of the EORs — are due to 
differential radial drift from a single reservoir, as there are sys- 



tematic chondrule size differences between the groups ( Clayton 



2003 Zanda et q/.|20 06), again consistent with S > 1 for them. 



7. Implications 

7.1. The value of S in steady disks 

To get a sense of the value of S and its variations in disks, we 
focus here on the stationary disk model (see Section|3]l. 

In order to estimate S , we will need a temperature prescrip- 
tion. Using equations (j7]i and ( A. 5 I, we have: 

.1/5 



T = max 



M 2 Q. 3 1 , 280K/ r /? 



1287T 2 cr SB k B (a} p 



-1/2 
AU 



,(36) 



with k the opacity, /?au = R/(l AU) and fi a dimensionless 
constant. The irradiation temperature is after Hayashi| (|1981|) 



but more realistic irradiation models (e.g. Chiang and Goldreich 
1997 ) should not yield too different power law exponents. 
We have: 



S = 



St 



3jt p s acr s 
~2~ MQ 



(37) 



Importantly, explicit reference to a has vanished — even in the 
viscous heating-dominated regime, the dependence of the tem- 
perature on (a)p is weak, see equation ( 36 1 — and is replaced by 



the mass accretion rate, which is an astronomical observable. 
This means that S can be evaluated without knowledge of the 
level of turbulence and its spatial variations in the disk. From 
equation ( 37 1, one sees that S is an increasing function of R and 
of time (that is, decreasing M), as would have been the case in a 
constant a disk (excluding the viscous expansion phase), where 



7 The spread AR is bounded by this value not only because of time, but also 
because K'driftl generally increases with R, as noted by Youdin and Shu (20021. 
By incurring a quicker drift of the outer wing of the distribution compared to 
its inner wing, this tends to squeeze it. The balance with diffusion, Drr/AR ~ 
ARd\v&m\/dR, indeed yields the value given in equation 1351. 
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Figure 4: Plot of the value of 5 for chondrule- sized bodies (p s a = 1 kg/m 2 ) in 
a steady disk, as a function of R\\], for three values of the mass accretion rate. 
Same parameters as Figure|2] To crudely mimick the drop of opacity as solids 
evaporate, we limit the viscous temperature to 1500 K, hence the first break 
in the curves (except for M = 10~ 9 M s ■ yr -1 ). S increases with heliocentric 
distance (and with decreasing mass accretion rate). The second break reflects a 
change in temperature regime from viscous heating- to irradiation-dominated. 



Sal 1 , and which we expect on that basis to be a fairly gen- 
eral behavior. Numerically, we have: 

0.19 f^f 3 



S =1 



,lkg/m 2 



10" 



\M-») \<Q'}p5cm 2 /g 



1/5 



0.12/ 7 



Ra\ 



(38) 



with M„ 8 = M/(10" 8 Mq ■ yr _1 ). This is plotted in Figure^ We 
note that as the steady disk solution likely overestimates S of 
real disks far from the Sun, this likely underestimates S there. 
The heliocentric distance where S = 1 (plotted in Figure ^ is 



R 



AU 



M_8 rnin 
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1 kg/m 
p s a 



2\5/3 



{a) p 5 cm 2 /g 



10" 



1/3 
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(39) 



7.2. Chondrite diversity: Space or time? 

Given our above expectation that S should be an increas- 
ing function of heliocentric distance and time, our conjecture 
that CCs formed with S < 1 while EORs formed with S > 1 
leaves two possiblities: either CCs formed closer to the Sun 
than EORs at some given epoch, or CC formation predated that 
of EORs, or some intermediate solution. We note that in both 
cases, the fact that EORs are more 16 0-poor than CCs would 
be qualitatively consistent with models of O isotopic evolu- 
tion involving isotopic exchange with 16 0-poor water from the 
outer disk (e.g. CO self-shielding, see Young et a/.|200 8): such 
evolution would indeed have begun in the outer disk and then 
spread to the inner dislrl 



However, samples returned from comet Wild 2 by the Stardust mission 
have oxygen isotopic signatures similar to CCs (Ogliore et al. 201 1 1, although 



Figure 5: Plot of location of the 5 = 1 line as a function of p s a (recall that 
chondrale-sized bodies have p s a = 1 kg/m 2 ). Same parameters as Figure [4] 
As discussed in Section |67T| this also corresponds to the maximum range of 
(significant) outward transport and decreases as expected with increasing size. 



The spatial alternative is unlikely, for the following reasons: 
The heliocentric distribution of asteroid taxonomic classes sug- 
gests that EOR parent bodies (presumably S-type asteroids and 
related classes) dominate the inner part of the main belt, while 
those of CCs (C-type asteroids and related classes) dominate 
its outer part (e.g. Bur bine et a/.||2008] l. Also, many CCs dis- 
play evidence of aqueous alteration and thus original presence 
of ice — stable far enough from the Sun — on their parent body, 
while evidence for water in EORs is comparatively marginal. 
This runs counter to the spatial ordering that the conjecture on 
S would suggest if time was not a significant factor. 

We thus suggest that CCs accreted earlier than EORs on av- 
erage, as also advocated by |Cuzzi et al. (2003 1 and Chambers 
( 2006| l. This leads to the scenario sketched in Figure [6] It is 
worth pointing out that this chronological interpretation could 
also account for the spatial zoning of the main belt via an accre- 
tion efficiency effect. Indeed, the planetesimal formation rate 



may be taken to scale roughly like £Q (e.g. Weidenschilling 
2004). For evolved phases of the disk where the surface den- 
sity is low, this could have had a nonnegligible effect — bearing 
in mind that meteorite parent bodies probably formed over a 
wider range of heliocentric distances than the present-day main 
belt (e.g. Wal sh et a/.|2011 1 — , accounting for preferential ac- 
cretion of EOR parent bodies in the inner solar system. 

Our chronology is certainly consistent with the greater abun- 
dance of CAIs and amoeboid olivine aggregates (contempora- 
neous with the former, e.g. Sugiura and Krot 2007) in CCs than 
in EORs ( |Cuzzi et aZ.|2003||Chambers|2006| l, since these are the 
oldest solar system solids. As to chondrules, age ranges seem 
as yet similar in CCs and EORs (1-3 Myr after CAI formation; 
|Villeneuve et a l. 2009). This would not exclude, however, that 
they were accreted earlier in CCs than in EORs. 



comets are thought to have accreted in the outer solar system. This would be 
compatible only in case of an early accretion or a reduced isotopic exchange 
with the purported ls O-poor water in the outer disk compared to the inner disk. 
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Figure 6: Cartoon of the proposed scenario. We show four schematic epochs 
(which may actually partly overlap): an early epoch where CAI are formed 
and transported (epoch "I"; it would correspond to stages A-C of Fig. 1 of 
jjacquet et a/.|201 T) , the epoch of chondrule formation (epoch "II"), the epochs 
of accretion of carbonaceous chondrites (epoch "III") and non-carbonaceous 
("EOR") chondrites (epoch "IV"). We shade in blue those regions where S < 1 
and in pink those where S > 1 . A dashed line marks the position where S = 1 
and drifts inward as the disk evolves. Millimeter-sized solids produced in the 
5 < 1 region are not efficiently transported to the S > 1 region. 



At early times, 26 Al was a significant heat source in plan- 
etesimals and many of them likely underwent differentiation 
(Sanders and Taylor 2005). Inasmuch as materials precursor 
to many differentiated meteorites should thus have accreted as 
S < 1, they would be expected to exhibit affinities to known 
CCs. Oxygen isotopes and depletion in volatile elements (e.g. 
|Wanke and D reibus 1986) would be consistent with such a 
parentage (Meibom and Clark 1999) but other stable isotope 
data exclude a complete identification (Warren 201 1 1. The CCs 
we observe today could have escaped extensive thermal pro- 
cessing e.g. because of the original presence of water (Ghosh 
et al. 12006] ), or as part of the unmelted periphery of their parent 



body as suggested by paleomagnetism measurements (Weiss 
\et al.\200% 

7.3. Coexistence of CAIs and chondrules 

We now investigate possible genetic relationships between 
CAIs and chondrules/matrix in carbonaceous chondrites. 

We first emphasize that the contribution of the "CAI factory" 
to CCs cannot be reduced to the CAIs themselves, for the fol- 
lowing two reasons: (i) There is no likely conceivable barrier 
for the condensation of most of the major rock-forming ele- 
ments after CAI formation. Most condensable matter initially 
in the "CAI factory" must hence have eventually condensed, 
e.g. heterogeneously on the first condensates, (ii) There is 
no dynamical reason that CAIs should be transported outward 
preferentially to the remainder of the CAI factory matter. Ac- 
tually, gas drag would rather lead to the opposite effect, but if 
X « 1 there, as we propose, there should be limited decoher- 
ence between CAIs and the remaining matter. 

It is noteworthy that then, incomplete condensation would 
not be ascribable to "gas loss" in the form originally hypoth- 
esized by Wasson and Chou ( 1974[ ) because gas would have 



remained coherent with the solids. This leaves kinetic effects 
as a plausible alternative explanation. 

From simple mass balance considerations (Grossman 2010), 
CAIs likely constituted no more than ~8 % of the condensed 
matter from the CAI factory (assuming a water- and organics- 
free CI composition). This means that the few percent of CAIs 
reported in CCs imply that the "CAI factory" contribution was 
likely a few tens of percent! Thus CAIs should not be thought 
of as a population foreign to the host chondrite and introduced 
only at the time of accretion. The fact that the Al/Si ratio in 
CCs, while supersolar, becomes subsolar if the CAIs are (men- 
tally) subtracted (Hezel et al. 2008| ), corroborates that there 



was some coherence between CAIs and the other components. 
This would imply that CAIs coexisted with chondrule precur- 
sors, then chondrules, and some of them must have witnessed 
chondrule-forming events, as supported e.g. by the occurrence 
of chondrules with relict CAIs ( |Krot et al.\2009) . 

That CAIs and the surrounding matter remained coherent is 
not at variance with the fact that their environment evolved to 
match that in which chondrules formed. Initial isotopic hetero- 
geneities would have been averaged out and mixing with mat- 
ter from the outer disk or infalling on the disk would have di- 
luted systematic isotopic anomalies. Fo r example, the excesses 
in 54 Cr and 46 50 Ti measured in CAIs ( Trinquier et al. 2009 1 



could have undergone dilution with isotopically differentiated 
meteorite-like matter to lead to the observed whole-rock com- 
position of CCs. This is in the range suggested by CAI abun- 
dances. 



8. Summary and conclusions 

We have considered a standard turbulent disk. The turbu- 
lence parameter a was assumed not to vary by more than a 
factor of several in the vertical direction, except perhaps for 
surficial active layers (in presence of a dead zone), and Schmidt 
numbers were assumed to be of order unity (although we kept 
them explicit throughout). We have ignored possible large-scale 
nonaxisymmetric structures and the influence of planets. 

In this class of disk models, we have derived the continu- 
ity equation for a population of grains in the tight coupling 
limit (that is, a dimensionless stopping time St <k 1), in or- 
der to study their dynamics until agglomeration in > 1 cm ob- 
jects. We have only considered transport of grains within the 
disk and ignored the effects of feedback of solids on the gas 
and that of accretion/fragmentation processes on the disk-scale 
transport of millimeter-sized bodies. A possible accretion bias 
against specific chondrite components was deemed irrelevant to 
explain the composition spectrum of chondrites although this 
issue would deserve more investigations. 

The dynamics of solid particles subject to gas drag are found 
to be essentially controlled by the "gas-solid decoupling param- 
eter" S = St/a. For X « 1, particles essentially follow gas 
motions and are well-mixed vertically, while for S » 1, they 
settle to the midplane and drift inward faster than the gas. We 
find that a one-dimensional model is as yet most appropriate to 
evaluate net radial transport. 
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We have attempted to estimate S for the millimeter-sized 
chondrite components at the time of accretion for the different 
chondrites. We have conjectured that S < 1 for carbonaceous 
chondrites (CCs) and S > 1 for non-carbonaceous chondrites 
(EORs), and have put forward four arguments: 

(i) Since particles produced close to the Sun cannot be ef- 
ficiently transported beyond the 5 = 1 line, the presence of 
calcium-aluminum-rich inclusions (CAIs) in CCs argues for 
S < 1 for them. 

(ii) The large chondrule-to-matrix ratio in EORs is consistent 
with preferential concentration of chondrules at the midplane as 
a result of settling, which requires S > 1 for them. 

(iii) Evidence for coherence between components of different 
grain size in CCs suggests S < 1, whereas Mg/Si fractionation 
in EORs relative to solar abundances implies S > 1 . 

(iv) The smaller range of oxygen isotopic compositions of 
chondrules in individual EORs compared to CCs could result 
from limited radial mixing in a S > 1 environment. 

S has been evaluated for steady disks and is found to be an in- 
creasing function of time and heliocentric distance. The above 
conjecture would thus favor earlier accretion of CCs compared 
to EORs, for which we listed some cosmochemical and planeto- 
logical supporting evidence. Another implication is that CAIs 
belong to the same reservoir as part of the chondrule precur- 
sors for CCs but significant mixing between different reservoirs 
must have occurred before chondrule formation. 

In this study, we have not attempted to introduce new 
physics, and we have recalled above some possibly relevant 
physical ingredients that we have ignored. Instead, we have 
tried to base our reasonings on as generic a disk as possible, in 
view of the existing literature on the transport of solids. The 
point of the present paper is that the simplest disk and radial 
transport models appear to lead to a robust conclusion about 
the space-time ordering of chondrite accretion. Our treatment 
of chondrites has been similarly generic: In being content with 
a first-order division of chondrites in two super-clans, irrespec- 
tive of the specificities of the groups that they encompass, we 
have dealt in all rigor with schematic end-members. In par- 
ticular, there is no reason to doubt that the transition from the 
S < 1 to the S > 1 regime was gradual and possibly some 
chondrite groups, e.g. those showing mixed affinities to CCs 
and EORs in terms of the criteria of Section 2 (e.g. CR chon- 
drites, Kakangari-like chondrites...) might be representative of 
accretion at S ~ I. 
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Appendix A. Vertical thermal structure of an optically 
thick disk region due to viscous dissipation 

In the optically thick limit (for Z > 1 g/cm 2 ) and assuming 
local dissipation of turbulence, the vertical temperature profile 
is governed by the transfer equation ( Terque m|2008| l 



(A.l) 



for timescales longer than the thermal timescale (Q.a) 1 . If we 
assume that k is vertically constant, this can be rewritten as 



d 2 T 4 



21Q.K 
16cr SB 



(A.2) 



with X'(z) = f Q p(z')dz', which, assuming that ac] is vertically 
constant, may be integrated as: 



T = T n 



1 - 



1/4 



(A.3) 



with r m id the midplane temperature and T m the "surface" tem- 
peratur^] both being related by: 



27 



128cr SB 

The solution may be approximated as: 



Tmid = max 



\ 128mcr SB 



1/3 



(A.4) 



(A.5) 



Equation (A.3i justifies the isothermal approximation over 
most of the column. Even if T- m <K T^, the temperature is 
within 20 % of T m \& for 80 % of the column. Had we given a 
temperature dependence to k, we would not have obtained a sig- 
nificantly different T m \&. In the irradiation-dominated regime, 
isothermality is also satisfied in the disk interior (between 
the "superheated" dust layers of Chi ang and Goldreich|1997] l. 
While we have assumed that ac 2 s was vertically constant, vari- 
ations of a factor of a few (notwithstanding the possibility of 
layered accretion) would also have little effect on T m \d because 
of the high power with which it appears in the equation. 
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